Feasibility of a single-parameter description of equilibrium viscous liquid dynamics 
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Molecular dynamics results for the dynamic Prigogine-Defay ratio are presented for two glass- 
forming liquids, thus evaluating the experimentally relevant quantity for testing whether metastable- 
equilibrium liquid dynamics to a good approximation are described by a single parameter. For the 
Kob- Andersen binary Lennard- Jones mixture as well as for an asymmetric dumbbell model liquid a 
single-parameter description works quite well. This is confirmed by time-domain results where it is 
found that energy and pressure fluctuations are strongly correlated on the alpha-time scale in the 
NVT ensemble; in the NpT ensemble energy and volume fluctuations similarly correlate strongly. 
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The physics of viscous liquids approaching the glass 
transition continue to attract attention Basic prop- 
erties like the origins of non-exponential relaxations or 
of the non-Arrhenius viscosities are still controversial. A 
question that is not currently actively debated is whether 
a single "order" parameter is enough to describe glass- 
forming liquids and the glass transition Q]- For more 
than 30 years the consensus has been that with few excep- 
tions more than one parameter is required, a conclusion 
that scarcely appears surprising given the complexity of 
glass- forming liquids. 

The prevailing paradigm in glass science regarding_''or- 
der" parameters may be summarized as follows [ll,H, H, 
S, H, 0, S Si ■ If is the drop in isobaric specific heat 
per volume going from liquid to glass, Akt and Aoip 
the same changes in isothermal compressibility and iso- 
baric thermal expansion coefficient respectively, and Tg 
the glass transition temperature, the Prigogine-Defay ra- 
tio n is defined by 
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Davies and Jones in 1952 proved that if there is just one 
"order" parameter, then 11 = 1 4). In their formulation, 
if a liquid is described by a single parameter, its lin- 
ear relaxations are all simple exponentials 4]. The vast 
majority of reported Prigogine-Defay ratios obey 11 > 1 
PI, which following Davies and Jones is consistent with 
the observation that linear relaxations are virtually never 
exponential. The case for requiring more than one pa- 
rameter got further support from the classical cross-over 
experiment of Kovacs [lO| as well as from several other 
experiments, all showing that glass structure cannot be 
completely characterized by a single parameter 11 , l2| • 
Part of this conventional wisdom may nevertheless be 
challenged. First, note that 11 is not strictly well de- 
fined. This is because to evaluate Acp, etc., one must 
extrapolate measurements in the glass phase into the liq- 
uid region, but the glass phase is not unique and it re- 
laxes with time. Secondly, as shown by Goldstein and by 



Moynihan and Lesikar (but perhaps not generally 

appreciated), it is possible to have systems with 11 = 1 
described by a single parameter with non-exponential dy- 
namics. Finally note that, although it has been conclu- 
sively demonstrated that not all aspects of glass structure 
can be described by a single number, this does not rule 
out the possibility that a single parameter is sufficient 
for describing a more limited range of phenomena, e.g., 
the linear thermoviscoelastic properties of a glass- forming 
liquid 14 , 1^ 1^ . It is this possibility we inquire into be- 
low, where all simulations were performed in metastable 
equilibrium with no reference to the glassy phase. 

Experiments carried out the last few years by Richert 
and Weinste in I l7l| . by Ngai, Casalini, Capaccioli, Paluch, 
and Roland [18[, and by our group [19|, do indicate that 
a single-parameter description may be appropriate for 
some situations. The direct motivation of the present 
work is a paper from 2004 by Mossa and Sciortino who 
simulated aging of a molecular model of ortho-terphenyl 
[20| . For small temperature steps they found that the 
location of the aging system in configuration space can 
be traced back to equilibrium states, implying that (20j 
"a thermodynamic description based on one additional 
parameter can be provided" for cases of non-linear relax- 
ations fairly close to equilibrium. This suggests studying 
the metastable equilibrium viscous liquid itself in order 
to investigate whether - and how - "one-parameter-ness" 
may be reflected in the equilibrium fluctuations. 

A single-parameter description of the metastable- 
equilibrium liquid dynamics by definition applies in the 
following situation [16l|. Suppose a linear-response exper- 
iment is performed where "input" infinitesimal tempera- 
ture and pressure variations ST(t) and 6p{t) are imposed 
on the liquid and "output" infinitesimal volume and en- 
tropy variations SV{t) and SS{t) are observed. Since 
highly viscous liquids exhibit time-scale separation be- 
tween the fast (vibrational) degrees of freedom and the 
much slower (relaxational) degrees of freedom, one ex- 
pects that each of the outputs may be written as instanta- 
neous couplings to the inputs plus a relaxing contribution 
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- the latter being a standard linear convolution integral 
involving the previoius temperature and pressure histo- 
ries. By definition, if the two relaxing contributions are 
proportional, the liquid is described by a single param- 
eter [l6j . In other words, a single-parameter description 
applies whenever a variable 5e exists such that 

5S{t) = 5e{t) ^ J^5T{t) ~ J^5p{t) 

5V{t) = 7 5e{t) + JTi5T{t) - J^Sp{t) . (2) 

In Ref. it was proved that in any stochastic dy- 
namics a single-parameter description applies if and only 
if the following "dynamic Prigogine-Defay ratio" is unity: 



Atp(w) 



(3) 



where the double primes denote the imaginary part of the 
response function. Moreover, it was proved that if the 
dynamic Prigogine-Defay ratio is unity at one frequency, 
it is unity at all frequencies. The relevant frequency- 
dependent thermoviscoelastic response functions are very 
difficult to measure - in fact no reliable measurements 
appear yet to exist. This single- frequency criterion could 
become useful for the interpretation of experiment, since 
it is likely that the first measurements will cover only a 
limited frequency range. In real experiments one can of 
course never prove that a number is exactly unity, but it 
seems reasonable to assume that the closer the dynamic 
Prigogine-Defay ratio is to unity, the more accurate a 
single-parameter description is. 

There is no simple physical interpretation of the pa- 
rameter Se except that it controls both entropy and vol- 
ume relaxations. In particular, there is no reason to 
believe that Se gives a complete characterization of the 
molecular structure in the way that the "order" param- 
eters of traditional glass science do. Despite this se- 
vere limitation it is not obvious that glass-forming liq- 
uids exist that are described by a single parameter to 
a good approximation. Recent results studying ther- 
mal fluctuations of less-viscous liquids give rise to op- 
timism, though: For a number of model liquids 21 1 - 
including the standard Lennard- Jones system - the po- 
tential energy and the virial correlate better than 90% in 
their thermal equilibrium fluctuations. Recall that the 
virial (when divided by volume) is the contribution to 
pressure from the molecular interactions, i.e., in addi- 
tion to the ever-present ideal-gas contribution deriving 
from the kinetic degrees of freedom [NUbT /V). The ob- 
served strong virial/potential energy correlations do not 
appear to depend significantly on viscosity. For highly 
viscous liquids, because of the time scale separation one 
expects that the slow (relaxing) contributions to pressure 
are given by the virial and that the slow contributions 
to energy come from the potential energy. Thus these 



variables should be highly correlated in their slow equi- 
librium fluctuations for systems similar to those studied 
in Ref. [ilj . This suggests that to a good approximation 
a single, relaxing parameter controls both quantities, in 
which case it is obvious to expect that Eq. ^ may apply 
as well. 

To test this the dynamic Prigogine-Defay ratio was 
evaluated for two viscous liquid model systems by com- 
puter simulation. The test requires that three frequency- 
dependent thermoviscoelastic response functions must be 
evaluated. For simulations carried out at constant vol- 
ume and temperature (the NVT ensemble) the relevant 
response functions are the isochoric specific heat per unit 
volume c„, the isothermal compressibility kt, and the 
"isochoric pressure coefficient" /3„ = {dp/dT)v- These 
quantities enter the dynamic Prigogine-Defay ratio [lij 
via 



ATy(w) 



(4) 



The fiuctuation-dissipation (FD) theorem [H [H, [H, [25 1 
implies that the above linear response functions may be 
determined from the equilibrium fluctuations of energy 
{E) and pressure (p) as follows. If AE{t) = E{t) — 
{E)v,T, Ap(t) =p{t) - (j>)v,T, and C^{f{t)] denotes the 
cosine transform of f{t) at frequency w, according to the 
FD theorem 



ATy(w) 



C^{(A^(0)A^(t))v,r}C^{(Ap(0)Ap(t))y,T} 

{C^{{AE{Q)Ap{t))v,T}f 

(5) 

If energy and pressure fluctuations correlate perfectly 
one has Ai?(i) cx Ap(i) and consequently ATy(w) = 1 
at all frequencies. Equation (O suggests that the dy- 
namic Prigogine-Defay ratio tests how well the two rel- 
evant quantities (here energy and pressure) correlate on 
the time scale defined by the frequency, it is sort of a 
time-scale dependent correlation coefficient squared. 

Figure 1(a) shows the frequency dependence of the 
imaginary parts of the three linear-response functions for 
the standard Kob- Andersen binary Lennard- Jones 80-20 
mixture [26| . The three response functions are very sim- 
ilar and all exhibit the characteristic asymmetry towards 
higher frequencies observed for real glass- forming liquids. 
Figure 1(b) shows the dynamic Prigogine-Defay ratio; as 
mentioned this number is unity if and only if the three 
imaginary parts are strictly proportional [l6l |. This is 
not the case, but the ratio is fairly close to unity in the 
frequency range of the main (alpha) process. Thus in 
the range of frequencies one decade above and below the 
loss peak frequency - the "alpha-relaxation range" - the 
ratio stays below 1.2. In accordance with the above inter- 
pretation of Ktv{'^) this indicates that for times in the 
alpha-relaxation range energy and pressure has a corre- 
lation coefficient larger than 0.9. 
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FIG. 1: (a) show the imaginary parts (scaled to maximum 
value) of —Cv{uj), — /3i,(aj) and l/Kriuj) for the Kob- Andersen 
binary Lennard Jones 80-20 mixture [2q|: (c) shows the same 
for a system of asymmetric dumbbells [271 ]. The dashed lines 
indicate Debye relaxation with relaxation times = 2000 
Can m/eaa aud Ta = 1600 ps respectively, (b) and (d) show 
the corresponding dynamic Prigogine-Defay ratios, Eq. (2). 
For both systems the total simulation time covers more than 
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FIG. 2: (a) and (b): Fluctuations of energy and pressure 
in the NVT ensemble for the Kob- Andersen binary Lennard- 
Jones mixture [1^ and the asymmetric dumbbell system [23| . 
Each point represents an average over a time interval of O.Itc 
where is defined from the loss peak frequency (Figs. 1(a) 
and 1(c)). Energy and pressure fluctuations are highly corre- 
lated, showing that a single-order-parameter description is a 
good approximation, (c) and (d): Fluctuations of energy and 
volume in the NpT-ensemble at the same state points as (a) 
respectively (b). 



We also simulated a highly viscous single-component 
molecular liquid. This is a system of "asymmetric dumb- 
bell" molecules defined as two Lennard- Jones spheres 
of different radii held together by a rigid bond 
As shown on Figs. 1(c) and 1(d) the imaginary parts 
of the response functions are similar, and the dynamic 
Prigogine-Defay ratio stays below 1.06 in the alpha- 
relaxation region. 

The results in Fig. 1 show that a class of viscous liquids 
exists where the dynamic Prigogine-Defay ratio is close 
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to one. This is consistent with earher computer simu- 
lations of the (poorly defined) standard Prigogine-Defay 
ratio n (Eq. ([1])) for various systems, showing that this 
quantity is often close to unity |28l. [29| . In Figs. 2(a) and 
(b) we plot energy and pressure equilibrium fluctuations 
in time for both systems. In order to focus on fluctua- 
tions in the alpha time range, both pressure and energy 
were averaged over one tenth of (defined as the in- 
verse loss-peak frequency^, corresponding to focussing on 
the inherent dynamics [30| appropriate for understanding 
the viscous behavior. The correlations are striking. Both 
systems were also simulated at constant temperature and 
pressure. Here energy and volume show similarly strong 
correlations (Figs. 2(c) and (d)). Thus as expected, the 
appropriateness of a single-parameter description is not 
ensemble dependent. 

To summarize, it has been shown that it is possible 
to investigate the single-parameter question by moni- 
toring thermal equilibrium fluctuations. The dynamic 
Prigogine-Defay ratio provides a convenient test quan- 
tity for this, a quantity that is also relevant for experi- 
ment because it refers to one single frequency. Clearly, 
several questions are still unanswered. For instance, it 
is not known whether there is any link to fragility; one 
traditionally expects that strong liquids to be single- 
parameter systems and fragile liquids to require a multi- 
parameter description. The findings of Rcf. [21,] . how- 
ever, indicate that network liquids are not well described 
by a single parameter in the above sense, whereas purely 
van der Waals liquids are. Clearly, more work is needed 
to clarify this and other issues, and we hope that the 
present results serve to encourage the physics and glass 
communities to reconsid the old question: "One or more 
"order" parameters?" 

This work was supported by the Danish National Re- 
search Foundation's Centre for Viscous Liquid Dynamics 
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